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1. Introduction 

The pressure of hot QCD, pqcd, constituting one of the most fundamental thermodynamic 
observables, has been under theoretical study for several decades now. Being of fundamental im- 
portance to cosmology (due to its influence on the cooling rate of the early universe) as well as 
of potential relevance to heavy ion collisions (through its influence on the thermodynamic expan- 
sion rate), pqcd has been computed using a variety of methods, including lattice Monte-Carlo, 
weak-coupling and large-Aff-methods, to name a few. 

Although at large temperatures T asymptotic freedom leads to the expectation that weak- 
coupling methods are sufficient to accurately describe the deconfined phase, it is well known that 
the infrared (IR) sector of QCD produces a challenge for perturbation theory, which for the case of 
the pressure arises at order g 6 T 4 [p. 

It has however been realized that this challenge can be overcome in an effective theory setup, 
where the problematic sector is described by a dimensionally reduced theory [^]. The key ob- 
servation is that thermal (equilibrium) QCD possesses three distinct physical scales, two of them 
generated dynamically. The contributions to the pressure (and to any other thermodynamic observ- 
able) from each of these scales can be obtained from carefully constructing and matching a series 
of effective theories (for details, see e.g. [^, |]]). 

The theories under consideration are (hard) QCD, electrostatic QCD (EQCD) and magneto- 
static QCD (MQCD), governing physics on length scales l/T (the typical scale for a particle in 
the heat-bath), 1 /gT (the dynamically generated screening length for longitudinal gluonic excita- 
tions) and l/g ' T (the dynamically generated screening length for transverse gluonic excitations), 
respectively. While the first two are amenable to perturbative calculations, MQCD is purely non- 
perturbative and has to be treated on the lattice. Viewing the gauge coupling g(T) as parametrically 
small (which is certainly justified at asymptotically high temperatures), these three scales are well 
separated, and can hence be dealt with individually via the effective theory setup. Schematically, 
for the pressure one can write pqcd = Pe + Pm + Pg> where each contribution depends on the 
matching scales. This scale-dependence will cancel in the sum, rendering pqcd a well-defined 
physical observable. 

In this letter, we will give a somewhat condensed account of what we currently know about 
the different contributions to pqcd- 

2. Status of the QCD pressure 

Below, we specify the contributions to the MS pressure pqcd = Pg + Pm + Pe from each 
physical scale individually, for the case of gauge group SU(/V C ) and Nf quark flavors. We will work 
at zero quark masses m q . = and vanishing chemical potentials \if = 0, and display all dependence 
on the MS scale p. 2 = A7ie~ y °}JL 2 by L = ln^r. Effects due to finite quark masses [||, ^] and 
chemical potentials ^ are available in the literature, but will not be discussed here. 

2.1 Contributions from the ultra-soft scale g 2 T, i.e. from MQCD 

Ultra-soft physics is not accessible by perturbative methods, due to the unscreened transverse 
gluonic sector, which would lead to severe infrared problems [[I]]. This sector is governed by a 
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three-dimensional pure gauge theory. Its only parameter is the dimensionful 3d gauge coupling g 2 M , 
which we write as = j^rf- The screening length gets generated non-perturbatively, making a 
numerical lattice Monte-Carlo treatment necessary. The detailed setup for how to incorporate the 
ultra-soft contribution into the physical pressure by a carefully defined mixture of perturbative and 
non-perturbative coefficients is explained in detail in Ref. The result is 



+ [pert] - [nspt] + [non-pert] + 0(e) 



(2.1) 



where oIa = N c z — 1, and OCg = §g — ^t^k 2 []4j, is a perturbative 4-loop coefficient. The three 
coefficients enclosed in square brackets originate from measuring the 3d YM pressure on the lattice 
and matching the result to the MS scheme. To be more precise, they are the following. 

• The first number stems from a non-perturbative lattice Monte-Carlo measurement of the 3d 
plaquette in pure SU(A^ C ) theory [^], 

[non-pert] = Hm jjS 4 ^1 - i-TrP^ - [ Clj 6 3 + c 2i 8 2 + c 3i 6 +c 4 ln/3] | (2.2) 

= 10.7(4) at/V c = 3 , 

where B = denotes the dimensionless lattice coupling, and ci 4 are divergences of the 3d 

r g M a 

lattice-regularized plaquette which can be computed in lattice perturbation theory. They read 
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The number in C2 is a sum of typical 2-loop (infinite-volume) lattice integrals [J10|, | 

2 (I? - 2 



3V4 >*-T +A * + h 



5.25449 



with (K is the complete elliptic integral of first kind) [ |11| , O , g] 
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1.013041(1) . 



(2.3) 

(2.4) 
(2.5) 

(2.6) 
(2.7) 

(2.8) 
(2.9) 
(2.10) 

(2.11) 
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The coefficient C3 has been estimated by numerical stochastic perturbation theory (NSPT) for 
N c = 3 [13] and, with higher numerical accuracy and full N c dependence, computed from 3-loop 



diagrams in lattice perturbation theory []14J]. In principle, it would be nice to know the full N c - 
dependence of Eq. fl2.2p. 

• The second number stems from an estimation of (the sum of all) 4-loop vacuum diagrams in 



lattice perturbation theory by NSPT [15], with the IR divergence regulated by massive gluon- and 



ghost-propagators (mass term \A and m cc in the action) [16]. NSPT works on a finite lattice of 



volume (aL) 3 , so the infinite-volume limit has to be taken first to ensure that the IR is regulated by 



the mass term only. The result is [17, 18] 



[nspt] 



(4tt) 4 
Zd A N c 6 

(4tt) 4 
Sd A N c 6 
47T 4 
729 



lim lim 



— TrP ) 



!-4 



C4 In ■ 



term 



am 



(2.12) 



+ c' 42 N c 4 + c' 43 N c 



25.8(8) = 13.8(4) in Feynman gauge at N c = 3 



Note that the errorbar of this number matches the precision obtained for [non-pert] . It would be 
nice to know all four coefficients, in Feynman gauge, either by a direct diagrammatic evaluation, 
or by doing NSPT for (at least) four different values of N c . 

• The third number stems from a matching 4-loop computation in the (3 — 2e)d continuum 
theory, regulated in the IR by gluon- and ghost-masses, with gauge parameter t, . Gauge depen- 
dence, introduced by the IR regulator, is guaranteed to cancel against that in [nspt]. The result 
reads 

-/) |^ + [pert] + ^(£)} (2.13) 

where J is the 1-loop massive tadpole integral J \/{p 2 + m 2 ). We choose Feynman gauge E, = 1 
which here leads to modified propagators 1 / /? 2 — ^ 1 / (p +m ), and obtain [19] 



[pert] = -3.73134481146281478501 in Feynman gauge 



(2.14) 



where the number can be expressed in terms of 18 fully massive 4-loop scalar master integrals 
[ ^0| , pT| |. For general we would have to calculate vacuum diagrams with two mass scales (ra 2 
and <^m 2 ), which presently is beyond our computational capabilities. 

The matching condition for the 3d gauge coupling reads [22, 23,|24] 
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For the g 6 pressure, only the leading coefficient is relevant. 



2.2 Contributions from the soft scale gT, i.e. from EQCD 

Soft-scale physics is governed by a three-dimensional gauge theory, coupled to an adjoint 
Higgs field. This adjoint Higgs theory possesses a small number of dimensionful coupling con- 
stants, which are related to the parameters of full QCD (being g 2 and T) by the equations given 
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below. The contribution of this sector to the pressure is given by 



Pm(T) ? t 4 J ■ • 
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(2.16) 



with the 4-loop coefficients a M = f§ - n 2 , p M = - §g - f§ In 2 - f In 2 2 + ^ 7i 2 - §^ k 2 In 2 + 

C(3) + 710 = —1.391512 p5], where 710 is the leading coefficient of a finite 3d scalar 4-loop 
integral that is known numerically only [25], 
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and the matching parameters are [ 1 1 , ||] 
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where we have used the beta-function coefficients j8o = 1 1 3 2z , j8i = ^ — -yz — z« and, for brevity, 
set z = Nf/N c . Writing Z n = £'(-«)/£ (-«), the coefficients read [g, [26|, ^ 

2 + z 



a E 4 



a E5 = 2a E4 Zi + -(l-21n2) , 
o 

1 - Z 
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(2.23) 
(2.24) 

(2.25) 
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a E7 = 2j3o7o + ---zln2, 



(2.26) 



as well as [28, 24] 



1$ = 4 j 8oa E4 (27o+Z 1 ) + ^(20 + 29z + 2z 2 )-2z( ? i + 31n2)-^z 2 ln2 
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(2.27) 



(2.28) 
(2.29) 

(2.30) 



where the 7, are expansion coefficients of the Zeta function £(1 — e) = — i + £" = o 7« (note that 
7o = 7£ = 0.577216). For the g 6 pressure, the g 6 terms of Eq. ( [2.20D are irrelevant. 

2.3 Contributions from the hard scale 2nT, i.e. from thermal QCD 

Hard-scale physics can be treated perturbatively, in a simple g 2 -expansion, without the need for 
resummations, thermal masses, or hard thermal loops. This is due to all IR effects being properly 
incorporated into EQCD and MQCD, and is one of the main conceptual advantages of using the 
effective theory setup. The contribution to the pressure from hard momentum scales reads 
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with ideal-gas coefficient a E i = 1 + \j, «E2 = — f (4 + 5z) [||, and [|3lj] 
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(2.32) 



and unknown coefficients j3 E i, which can be determined e.g. by a 4-loop computation of vacuum 
diagrams in thermal QCD. Since pqcd is physical, the divergent and scale-dependent parts of j3 E i 
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are related to the other coefficients introduced in the above, serving as a valuable check on this 
open computation. Specifically, from 2-loop running of the 4d gauge coupling 



I 2 = = sHfr) +| 4 (A>)(-2A>^) +| 6 (Mo)(4j3 2 ^ 2 - 2fr£) 
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The remaining g 6 -coefficient, j8 E i however, entails a four-loop computation of all connected vac- 
uum diagrams involving quarks, gluons and ghosts, a computation that has so far not been tackled 
due to the formidable task of solving many genuine 4-loop sum-integrals. From diagrammatic 
arguments, it is clearly a polynomial in z = Nf/N c , 



=#o + z#i+z 2 #2 + z 3 # 3 , 



(2.37) 



and we will in the following indicate how two of its coefficients (the first and last) can be crudely 
estimated numerically already. 



3. Putting everything together 



Expanding in e, all poles cancel, as they should. In practice we make use of Eqs. ( |2.15| ),( |2.19| ) 
and ( 2.20| ) to re-expand all terms with a factor l/e or L in Eqs. (2_A) and ( 2.16 ) in terms of g 2 . 



After cancellation of the poles (and taking into account terms like | • e), we can now take the limit 
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Collecting explicit logarithms L, they precisely cancel the scale dependence of g 2 up to the 
order of the computation, and can hence be absorbed by writing 



g 2 = g z + t2p L + g b (4/^Z/ + 2j3iL) + 

Note that this coupling is explicitly scale independent to the order we are working, d ln ^2g 2 = €?(g & ). 
We now have the full pressure as a sum of its ultra-soft, soft and hard parts as 
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(3.4) 
(3-5) 
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Figure 1: Left panel: The normalized QCD pressure pqcd/psb at Nf = plotted versus the effective 
coupling g from Eq. (3.3). The g 6 coefficient depends on an unknown parameter A as defined in Eq. (3.8), 
and the different curves correspond to choosing A = —4000 (lowest curve) to A = +10000, in steps of 2000. 
Left panel: The same, plotted versus In The black dots correspond to lattice data from 



~3 1 , ~2 ~2 
m E ^+g E m E 



ln(2m £ ) - - 
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+ gE m E 
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Ph = &E1 +g &E2+g [&E3 ~ 180a E5 ] +f |0 m - 180 (0E2 + a E 4/3 E3 + 

The g 6 coefficient of /?qcd hence depends on a constant 

A = j3 E1 ±7205 NP ±3855 NS PT, 



(3.6) 
• (3-7) 

(3.8) 



where we recall that j3ei stands for the result of the open 4-loop computation, and the 8 parameter- 
ize the error-bars of the numerical constants from Eqs. ( 13|J2.13 ) as [10.7 ± 5np] and [25.8 ± 5nspt], 
respectively. From above, we have 5np = 0.4 and 5nspt = 0.8, so A = j8 E i ± 600. In the following, 
we will for simplicity set 5np = 5nspt = 0, remembering the induced error-bar on /3 E i. Using the 
same coupling as in p h , the above matching conditions now read 
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gE 
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r«E4+| 4 aE6 + 



f + g 4 a E 7+g 6 YEi + 
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(3.9) 
(3.10) 
(3-11) 
(3.12) 



We would now like to plot the result for pqcd- Identifying the non-interacting (ideal-gas; 
Stefan-Boltzmann) limit as psb = <^4T 4 |^G! E i, we could display the normalized pressure /?qcd//>sb 
as a function of the coupling g, for fixed Nf. This is done in the left panel of Fig. |J, at N{ = and for 
various A. Our goal, however, should be to try to make contact to existing lattice determinations of 
the full pressure, where typically /jqcd / Psb is given as a function of T /T c . Continuum-extrapolated 
lattice data exist for A*f = only, so in the following we will restrict to this special case. Aiming 
for this rather phenomenological comparison, we evidently need to make some choices, specified 
below. 
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Figure 2: Left panel: The two real branches of W (z) versus z. The upper (dashed) branch is Wq(z), the lower 
(solid) branch is W_i(z). Right panel: The effective coupling g from Eq. (3.3) plotted versus In ^, using the 
choices explained in Sec. 3. The upper/lower curve corresponds to the uncertainty in scale choice stemming 
from fixing p. and determining T c /A^ with (5m, &r c ) — (2.0,0.9)/(0.8, 1.1), the bigger effect coming from 
the latter parameter. 



We use the running 4d coupling from the exact solution of the 2-loop RGE equation, 



a2 



(A) 



I +\V i [ -^exp 



.l_ 2 |(L + lng) 



(3.13) 



Here, W_ i (z) is one of the two real branches of the Lambert W function (see e.g. [33] and the left 
panel of Fig. |2[ W(z) is the function that satisfies Wexp(W) = z). Note that the above solution 
entails two choices: The branch of the W -function and the integration constant were chosen in 
accord with asymptotic freedom (note that the argument of W — > 0~ for p. — > °°) and the 'usual' 
definition of A^s (being the absence of a 1/ln 2 p, term in the asymptotic expansion of g(jX) at large 



H). Indeed, at large L 



In -A- the expansion W-i(— e) 

A MS 

Pi ln(2L) 



ln£ - lnln| + ^(1/lne) reproduces 



+ GiL^) , in accord with e.g. Ref. \M\. 



fm = l/(2M+|ln(2L) + ^(L-)) = ^ L 4i , 2 

Although in principle all dependence on the renormalization scale pL, entering through L, is 
of higher order, in practice we need to fix it once we need numerical values for the coupling g. 



Following [g4j], we choose the scale p. by the principle of minimal sensitivity applied to the 1-loop 
result for g\, and then estimate the scale-dependence by a variation of a factor of 8^ = [0.8..2.0] 
around this u opt , obtaining L = — ^f- -t-ln-jf-. The slightly asymmetric choice of 8 U here reflects 

the fact that the 1-loop g 2 E falls off more steeply on one side of the plateau than on the other. 



To compare with continuum-extrapolated lattice data []32|], we use = 1.225^ where 8t l . 



[0.9..1.1] encompasses the central values and error bars of estimates of this quantity from different 
lattice collaborations (for a summary of the different methods and results, see [|24|]). This would 
translate into a horizontal error bar for the lattice data when plotted against T/A^g. 

In the right panel of Fig. 0, we have plotted the effective coupling g as defined in Eq. (3.3), 



converted to a function of In Note that its value is smaller than 0.2 even at T c . 

The normalized pressure pqcd/psb, converted to a function of ln^- along the lines above, is 
displayed in the right panel of Fig. For comparison, the continuum-extrapolated lattice data of 
Ref. [32] has been included as black dots. The figure suggests a value for the Nf = coefficient of 
the unknown constant j8ei , #o ~ 6000 ± 600, bearing in mind the error bar defined in Eq. (3.8). 
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4. Discussion 

We have currently no idea what the 4-loop hard-scale coefficient j8ei is, even though it can be 
computed diagrammatically. As already mentioned above, it should be a polynomial in z = Nf/N c , 
/fei = #o + +z 2 #2 + z 3 #3, where only a single (ring) diagram contributes to #3, suggesting it as 
the first test-case for 4-loop sum-integral technology. 

It seems possible to give an estimate of the highest-Aff contribution to Pei from the large-Ab- 
solution for the pressure, since terms of order g 6 Nf 3 originate from the hard-scale pressure only. 
In [35] this was attempted by fitting the numerically known exact large-Aff pressure with a polyno- 
mial in g. This results in #3 = ^ [+20(2)] - f ln2 (l + 12y - ^ln2 + I6Z1 - 8Z3) » [+36(1)], 
where the terms proportional to In 2 originate from translating the choice of renormalization scale 
p, = 7iT of [ [351 ] to our definition of jSgi, where powers of L = In were subtracted out. 

Furthermore, fitting the full g 6 pressure at Nf = and N c = 3 to lattice data around 4T C [32] 
suggests a value #0 ~ 6000 ± 600, if one takes the conjecture for granted that all higher-order 
corrections sum up to a subdominant contribution. There is no guarantee whatsoever that this 
conjecture holds, making a perturbative computation of the #, unavoidable. We take the above 
check against the lattice data as indication that the effective theory setup has a chance to analytically 
describe the transition from temperatures as low as a few times T c to infinite temperatures, in terms 
of computable corrections to the ideal-gas limit. 
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